Heavy Quasiparticles in the Anderson Lattice Model 
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An exact-diagonalization technique on small clusters is used to study the dynamics of the one- 
dimensional symmetric Anderson lattice model. Our calculated excitation spectra reproduce key 
features expected for an infinite Kondo lattice such as nearly localized low-energy spin excitations 
and extended regions of 'heavy-quasiparticle' bands. We show that, in contrast to the hybridization 
picture, low-energy spin excitations of the nearly-localized /-electron system play a key role in the 
formation of an almost dispersionless low-energy band of heavy quasiparticles. 

PACS numbers: 71.27.-ha, 71.28.-hd, 75.20.Hr 
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The origin of the anomalous behavior of /-electron 
compounds is an unresolved issue in the theory of 
strongly-correlated electron systems. Thereby the way in 
which a periodic array of magnetic ions interacting with 
a sea of conduction electrons can give rise to either the 
extreme low-energy scale in the Landau-type quasiparti- 
cle bands of heavy-Fermion compounds [y^] or to gaps 
of apparent many-body origin in the excitation spectra 
of Kondo insulators is not yet understood. On a 
phenomenological level, heavy-Fermion compounds have 
been described with considerable success by the 'renor- 
malized band theory' [Q where the effect of electron cor- 
relations is described by the renormalization of on-site 
energy and hybridization strength of the magnetic ions. 

In this Letter we study the Anderson lattice model 
(ALM), the simplest model relevant for /-electron com- 
pounds, by Lanczos diagonalization of small clusters 
and show that such a renormalized band picture on one 
hand may provide a reasonable phenomenological de- 
scription of the dispersion relations, but on the other 
hand is not really adequate on a microscopic level: con- 
trary to the one-particle picture, the heavy quasiparti- 
cles may be viewed as loosely-bound states of conduction 
electrons and spin-wave-like excitations of the nearly- 
localized /-electron system. The emerging picture is thus 
more reminiscent of the spin polaron discussed recently 
on the basis of a scmiclassical treatment of the Kondo lat- 
tice . We would like to stress that due to the small size 
of the clusters our calculations can neither reproduce the 
exponentially small energy scales present in the Anderson 
impurity problem nor do they allow for the derivation of 
an accurate phase digram. Nevertheless one may expect 
that the relative magnitudes of energy scales and hence 
the nature of the low lying states are reproduced correctly 
by our calculations; in that sense our results for dynam- 



ical quantities are complementary to the scaling theories 
in Ref. Q or the renormalization group calculations for 
the ground states of larger systems in Ref. |^ . 

We consider a tight-binding version of the one- 
dimensional ALM defined by the Hamiltonian 



H= - t (4c,. + H.c.)-^^(cL/.. + H.c.) 

< ij > cr ia 
i 

where Cia- {fia) is the annihilation operator for an elec- 
tron of spin cr at site i in the c- (/-) orbitals and 
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fiafic- Model parameters are hopping strength t 
between nearest-neighbor c orbitals, mixing V between 
c and / orbitals, on-site repulsion U at the / orbitals. 
The on-site energy of the / orbitals is taken to be — C//2, 
i.e., we consider the 'symmetric' case. We focus on elec- 
tron densities close to 'half filling', i.e., 2Ns electrons 
in Ns unit-cells (a unit cell contains one c and one / 
orbital). Restrictions on memory space and computer 
time necessitate to choose Ns<Q. To get additional in- 
formation, we employ twisted boundary conditions [p|jTc|] 
(BC). We require CAr^+i^.^e^'^cio- and f^^+i^^e^'fifia by 
introducing an arbitrary phase ip; the allowed momenta 
are then k={2Trn+(p)/Ns with n=0, . . . , A'^s — 1. It has 
been pointed out 12| that in the half-filled case finite 
size effects can be minimized by choosing ip={Ns/2) ■ tt. 
As will be seen below, combination of spectra obtained 
with different values of tp gives remarkably smooth 'band 
structures'. However, none of our conclusions to be pre- 
sented below depends on such an assignment of bands. 

We first consider the single-particle spectral function 
A^{k,LLj) defined as a sum Aj{k,uj) = A~{k,—uj) + 
A^(k,uj) of the photoemission (PES) spectrum 
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and the inverse photoemission (IPES) spectrum 
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where (|^^)) denote the ground state energy (wave 
function) with N electrons and twisted BC of phase 
if. The operator 7^0- refers to the Fourier transform 
of the operator for either conduction electrons (cia) or 
the / electrons {fia)- Results for A^{k,uj) obtained by 
the standard Lanczos procedure are given in Figure 1. 
One can identify the 'upper- and lower-Hubbard bands' 
for the / electrons, separated by an energy ~J7. They 
are dispersionless and somewhat broadened, with almost 
pure / character. In addition to this typical strong- 
correlation feature, there are two bands which are more 
reminiscent of noninteracting electrons: one can iden- 
tify the unrenormalized c-electron band of width At, ap- 
parently split into two bands by hybridization with a 
'renornialized' / level in the middle of the Hubbard gap. 
This feature first results in a well-defined gap between 
PES and IPES spectra and second in extended regions 
of 'heavy' bands with apparently pure / character. The 
dispersion of the two 'hybridization bands' as well as the 
change from almost pure c character to almost pure / 
character around 7r/2 thereby are both roughly consis- 
tent with the picture of noninteracting electrons. The 
spectral weight of the parts with / character, however, 
is substantially smaller than in the parts with c charac- 
ter; comparison shows that this asymmetry is the more 
pronounced the larger the ratio U/V. 

The change ofAj{k, uj) with hole doping is at first sight 
completely consistent with the picture of noninteracting 
electrons [Q: the chemical potential seems to shift into 
the 'heavy' band, so that a kind of Fermi surface emerges, 
and upper- and lower-Hubbard bands remain unaffected. 
In addition to this rigid-band-like behavior, however, 
there is also a modification of the 'light' parts of the band 
structure, far from Ep: c-type spectral weight is trans- 
ferred from PES to IPES near 7r/2, i.e., the Fermi mo- 
mentum for a half-filled band of unhybridized conduction 
electrons. The change of Ac{k,uj) thus is reminiscent of 
unhybridized conduction electrons. On a phenomenlogi- 
cal level, this could be reproduced if one assumed that the 
'renormalized' /-level energy is pinned near the chemical 
potential of N—N^ unhybridized conduction electrons, 
i.e., the Fermi energy of a 'frozen-core' band structure. 

We now want to clarify the nature of the heavy-band 
states. Important information can be obtained from the 
momentum distribution function nj„{k)={'fl.^^kcr), more 
precisely, we study the change of n^^{k) upon remov- 
ing one electron. In a six unit-cell system, we evaluate 
the difference An-yo-(^) between the nj„(k) in the lowest 
state with 5 down-spin and 6 up-spin electrons at the to- 
tal momentum /ctot and the n-y„(k) of the ground state 



at half-filling. We choose /ctot such that the single-hole 
state belongs to the 'heavy' part of the band. In the hy- 
bridization model, the creation operators in the lower hy- 
bridization band would read = Ukcl^ + Vkfl^, so that 
AneT(fc) = 0, A7iex(fc) = -\uk\^Sk,-k,^,, An/T(fc) = 0, 
and Anfi{k) — —\vk\'^Sk,-ktof Since one may expect 
Uk~0 and Vk'^l in the heavy band, the electron is re- 
moved only from the / species with spin down and at 
k=—ktot- The calculated results for An^o-(fc) are shown 
in Figure 2, where we note the following features, almost 
all of which are in contrast to these predictions: 

(i) Independently of the actual momentum /ctot of the 
single-hole state, c electrons of both spin directions are 
removed at the two fcp, 

(ii) The resulting loss of up-spin electrons is compensated 
by an almost /c-independent spin polarization of the / 
electrons, 

(iii) As the only agreement with the hybridization model 
there is an extra 'dip' in nf^(k) for fc=— fctotj which how- 
ever diminishes rapidly in magnitude for decreasing V/U. 
These results establish first of all that the 'heavy quasi- 
particle' is predominantly a 'missing c electron' with only 
small admixture of / character (for large U/V). By con- 
trast, the pure / character of the lower Hubbard band 
suggests that it is in this band where an / electron is 
missing. We thus have an energy separation of c-like and 
/-like degrees of freedom (of order U/2), in contrast to 
the hybridization scenario. However, there must be some 
mechanism which renders the missing c electron 'invisi- 
ble' in Ac{k,uj) when k is in the 'heavy' band. 

As for this latter issue, we note that the spin polariza- 
tion of the / electrons suggests the presence of a spin ex- 
citation. We therefore consider the spin-excitation spec- 
trum 
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where is the Fourier transform of either the total-spin 
raising operator c\^Ci[ + fj-^ fii (a=tot) or the /-electron 
spin raising operator f}^fii (a—f). The calculated spec- 
tra (see Figure 3) show strong low-energy peaks with neg- 
ligible dispersion, which probably are the (almost) local 
singlet-triplet excitations expected for a Kondo lattice. 
In Sf{q, Lo), these low-energy peaks are enhanced whereas 
the smaller peaks at higher energies are suppressed: ob- 
viously the spin-flip of an / electron in the ground state 
to excellent approximation produces another eigenstate. 
There is a pronounced k dependence of the peak inten- 
sity, similar to spin waves in an antiferromagnet. One 
may assume that this reflects the antiferromagnetic spin 
correlations due to the RKKY-type interaction. We also 
study the change of the momentum distribution of the 
half-filled system due to a spin excitation; more precisely 
we consider the difference between the momentum distri- 
bution for the lowest state with 5^=1 and momentum tt 
(i.e., the final state for the low-energy peak in iS'q((J, w)) 
and that for the ground state. This difference is shown 
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in the inset of Figure 3. Whereas the c electrons remain 
virtually unaffected by the spin excitation, there is an 
almost /c-independent polarization of the / electrons, as 
one would expect it for a quantum spin system without 
charge degrees of freedom. Again we find a remarkable 
degree of separation of the c- and /-electron 'subsystems', 
which may also provide a natural explanation for the 
strongly difi'erent spin and charge excitations found in 
previous studies ,^2|jl|J|,|l|-^ of Kondo insulators. 

Let us now combine the above results to obtain a sim- 
ple picture of the heavy states. Since they represent the 
parts of the PES /IPES spectrum with the lowest excita- 
tion energy, let us consider the limit V—>0 and ask 'How 
can we remove or add an electron so as to lower the en- 
ergy most efficiently?'. In the half-filled ground state, 
there is on the average one /-electron per unit cell, with 
only a small admixture of the empty or doubly occupied 
/ site. Removing or adding an /-electron will on the av- 
erage raise the energy by U/2, and thus is unfavorable. 
Accordingly, the state fka\'^!^), which would be the most 
natural ansatz within the hybridization picture, has only 
small overlap with the true 'heavy state', particularly in 
the strong correlation case (i.e. small V/U). One mea- 
sure for the weight of this state in the ground state would 
be the 'depth' of the dip in nfc{k). On the other hand, 
a c-electron can be removed or added with practically 
no cost in energy if that is done near fcp. Next, the /- 
electron spin excitations with their small excitation en- 
ergies offer a way to dispose of 'excess momentum' with 
almost no cost in energy. This suggests to remove or 
add the c-electron always at k^, and transfer the excess 
momentum to an /-spin excitation. This picture imme- 
diately explains the reduction of 71/^(^5-)' ^ ^^^^ 
spin polarization of the /-electron system due to the ac- 
companying /-spin excitation. We are thus led to the 
following ansatz for a hole-like 'heavy state': 



-Cfce^5+(fc-f jl*^). (5) 



Here S'| [q] is the z-spin operator for the /-electrons with 
momentum transfer 5, and Ufe and Vk are (variational) pa- 
rameters. The state (Eq. (H)) has momentum k, z-spin 
1/2, and total spin S'=l/2, i.e., the spins of /-electron ex- 
citation and c-electron hole maximally compensate each 
other. This is reminiscent of the 'quenching' of a Kondo- 
impurity spin due to bound-state formation. The vari- 
ational parameters in (^ are determined from the re- 
quirement that \^{k)) has norm 1 and maximum overlap 
with the exact 'heavy' state with momentum k. Figure 4 
shows the overlap |(^'(fc)|^^^^)P for different values of 
V/t and U /t at fc=57r/6; here |^;^^^) denotes the exact 
'heavy' state and |^'(fc)) is given by Eq. (||). For compar- 
ision, the overlap of the state f-ki\'^^) (normalized to 
unity) with \^^^^) is also shown (in the hybridization 



picture, the latter quantity would be 1). While the 'bare 
/-electron' is a good approximation only in the small 
U /V case, the overlap of the state in Eq. (||) is >90%, 
for all parameter values, so that we find a good descrip- 
tion of the heavy-band states in the strong correlation 
region. 

In summary, we have studied the single-particle spec- 
tral function and dynamical spin-correlation function for 
finite clusters of the Anderson lattice model at and near 
half filling. Despite the necessarily rather coarse energy 
scales available in the clusters, our results do reproduce 
key features expected for infinite Kondo lattices, namely 
extended heavy bands and almost dispersionless low- 
energy spin excitations. On a phenomenological level, 
the low-energy parts of the spectral function can be de- 
scribed reasonably well by a renormalized band picture 
where an 'effective / level' pinned to the 'frozen-core' 
Fermi energy mixes with the conduction band. This pic- 
ture, however, has not much significance beyond a purely 
phenomenological level: there is a clear separation be- 
tween the low-energy 'hybridization bands' which cor- 
respond to a missing (or extra) c electron and the two 
'Hubbard bands' which correspond to a missing (or ex- 
tra) / electron. The heavy quasiparticles rather have 
the character of loosely bound states between a conduc- 
tion electron at the Fermi momentum of the unhybridized 
conduction-electron system and a spin-wave-like excita- 
tion of the /-electron lattice which acts very much like a 
pure quantum-spin system. 

Since the 'spin polaron bands' are formed by bound 
states, rather than by true hybridization, it seems natu- 
ral to assume that the breaking of the bound states will 
completely remove the heavy parts of the band struc- 
ture and leave behind only the frozen-core Fermi surface. 
Since the heavy quasiparticles involve the spin compen- 
sation of / excitation and c hole, it is moreover clear 
that they can be broken by a magnetic field. Then, the 
breaking of the heavy polarons by a magnetic field and 
the corresponding collapse of the Fermi surface to the 
frozen-core volume then appears as a natural explana- 
tion for the so-called metamagnetic transition associated 
with 'itinerant-to-localized' nature of / electrons |20[| . 
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FIG. 1. Single-particle excitation spectra Aj{k,Lj) for 
Ns=6 with different ip for V/t=l and U/t=Q. The spectra at 
fe=0, 7r/3, 27r/3, and tt correspond to (^=0 (i.e., periodic BC). 
The / spectra are multiplied by —1 for better distinction, the 
Lorentzian-broadening is e=0.02t. The upper panel shows the 
spectra for the half-filled ground state, the lower panel for the 
ground state with 5 up and 5 down electrons (i.e., with two 
holes). The vertical dashed line shows the chemical potential. 



FIG. 2. Difference An-yo-(fc) between zero-hole and one-hole 
ground states at U /t=Q; (a) V/t=l, fctot=57r/6, ip=-K (i.e., 
antiperiodic BC), (b) V/t=l, fctot=7r, ip=Q (c) l//t=0.5, 
A:tot=57r/6, ip=-n. As a reference n^a{k) for the zero-hole 
ground state at V/t=l is shown in (d). 

FIG. 3. Spin-excitation spectra Sc,{q,oj) for V/t=l and 
U/t—6, ifi=0. The Lorentzian-broadening e=0.02t. Inset 
shows the change of n~fcr{k) due to the spin excitation. 

FIG. 4. Overlap of the state given by Eq. (^ (white sym- 
bols) and the state |*^)/(*^|/l^^/_fei (black 
symbols) with the exact 'heavy' state l^;^"'^) at fc=57r/6 and 
ip=iT as functions of U /t. 
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